DIFFERENCE APPROXIMATION FOR LOCAL TIMES OF 
MULTIDIMENSIONAL DIFFUSIONS 
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Abstract. We consider sequences of additive functionals of difference approximations for uniformly 
non-degenerate multidimensional diffusions. The conditions are given, sufficient for such a sequence 
to converge weakly to a Vy-functional of the limiting process. The class of the ly-functionals, that can 
be obtained as the limiting ones, is completely described in the terms of the associated VF-measures 
^ by the condition 
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1. Introduction 

In the paper, we consider an M'"-valued diffusion process X defined by an SDE 

(1.1) X{t)=X{0)+[ a{X{s))ds+ [ b{X{s))dW{s), t G M+, 

Jo Jo 

and a sequence of processes Xn,n > 1, with their values at the time moments G N given by a 
difference relation 



A:-l\ k-l 



(1.2) Xn[-]=X„ 

' n J \ \ n J J n \ \ n 
and, at all the other time moments, defined in a piece-wise linear way: 
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Here and below, PV^ is a Wiener process valued in M™, {£,k} is a sequence of i.i.d. random vectors in 
M™, that belong to the domain of attraction of the normal law, are centered and have the identity for 
covariance matrix. Under standard assumptions about coefficients of the equations (jl.ip , ()1.2p (local 
Lipschitz condition and linear growth condition), the distributions of the processes Xn in C(M"^,R"^) 
with the given initial value X„(0) = x converge weakly to the distribution of the process X with 
X(0) = X (see P). 

In the paper, we deal with the following problem. Let {(/9'^'*,0 < s < t} be a M^-functional of 
the diffusion process defined by (II. ip . This, by definition, means that (see p]. Chapter 6) (/9 is a 
non-negative homogeneous additive functional with its characteristics 

{f\x) = ^[/'*|X(0) = x],t > 0,x G M™} 
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satisfying the condition sup^. /*(x) < +00, t £ We consider a sequence of non-negative additive 
functionals < s < i}, n > 1 of the processes X„, of the form 

(1-4) = (X„) i Yl ^"(^"(^))' 0^«<*' 

k:s<-<t 

— n 

and give the conditions sufficient for the joint distributions of conditioned by X„(0) = x 

(x S M"* is an arbitrary point), to converge weakly to the joint distribution of {ip,X), conditioned 
by X(0) = X. 

One motivation for posing such a problem is the following one. It is well known that the theory 
of additive functionals of M"*-valued Markov processes is closely related with the potential theory. 
There exists a one-to-one correspondence between VF-functionals and so called T^- measures (see j2], 
Chapter 8); every W^-functional ip can be written at the form 

(1-5) <^^'*= f 4^{X{r))dr, 0<s<t, 



where ^ is the corresponding VF-measure, A™' is the Lebesgue measure on W^. In general, W^- measure 
fj, is not absolutely continuous w.r.t. A™; for singular /i equality p.5|) is a formal notation, that can 
be substantiated via an approximative procedure with /i approximated by an absolutely continuous 
measures. The functional if, given by (jl.Sp . is called the local time of the process X, corresponding 
to the measure /i. 

Given ly-functional (p, one can construct the sub-process X"^ with its transition probability given 

by 

(1.6) PiX'^it) Gr|X'^(0) =x] =^[Ir(X(t)) •exp(-/'*)|X(0) =x], t E M+, x e M"', T G 5S(M™) 

(see [2j, §6 of Introduction); the measure /i is interpreted as the killing measure for this process. Let 
X be the diffusion process given by the equation (jl.ip . then its generator is equal 
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1=1 *ii=i 

The well known Feynman-Kac formula gives the generator of the process X'^ at the form 21'''/ = 
21/ — • /, i.e., for every continuous bounded function g{-), the solution to the following Cauchy 
problem for the second order parabolic PDE 



(1.7) u'tit, x) = ^ ai • -Q^u{t, a;) + 2 E "^^J ■ dx dx - ^^^' ~ Ix^^^^^^^' ^ 

i=l * iii=l * 

has the following probabilistic representation: 



(1.8) n(t, x) = E[g{X{t)) ■ exp{-ip'^^')\X{0) = x 



Let us note that the term in the equation p.7|) . that corresponds to the "heat flow-out", is a 
generalized function and equation (jl.7p should be interpreted in the generalized sense (therefore, it is 
natural to call this equation a singular one). On the other hand, its solution is a bounded measurable 
function due to the representation ()1.8p . 
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Put 

(1.9) Unit,x) = •exp(-(^0'*)|X„(0) = x]. 

The main result of the present paper (Theorem I2.ip provides convergence in distribution of the 
difference approximation {ipn,Xn,) to {ip,X); thus, under conditions of this theorem, 

Unit,x) ^ u{t,x), n ^ OO, tEM+,xeM'". 

This means that this theorem, in particular, gives an opportunity to apply, in a standard way, the 
Monte-Carlo method for numerical solution of the Cauchy problem for the singular parabolic equation 

Let us make a short overview of the bibliography devoted to weak convergence of the functionals of 
the type (jl.4p . For one-dimensional random walks and difference approximations of one-dimensional 
diffusions, there exists a large variety of limit theorems for the associated additive functionals. We 
do not discuss these results in details, since we are mostly interested in a multidimensional case, and 
refer to the monographs [S],!!] and papers [S] - [TU]. 

In the multidimensional case the situation is essentially different. The author does not know any 
paper where a limit theorem for the functionals of the type (II. 4p would be proved in the case where 
Xn approximate a non-additive diffusion process X (i.e, where the coefficients a, b are non-constant). 
The only multidimensional limit theorem, known for the author, for the additive functionals of the 
type (jl.4p . is given in the paper [11] in the situation where X„ is a multidimensional aperiodic lattice 
random walk and X is the Brownian motion in M"^ (also, the paper [12] deals with the closely related 
problems). The significant difference between the results available in the one- and multidimensional 
cases can be naturally explained by the fact that the structure of the class of H^-measures is much 
more complicated in the second case than in the first one. For the Brownian motion (and, also, for 
any non-degenerate diffusion), for m = 1, every finite measure is a VF-measure. For m > 1, any 
measure 6z,z € is not a VF-measure {4^ there does not exists the local time at any fixed point 4^ 
every one-point set has its capacity equal to zero 4^ every one-point set is a polar set). Therefore, in 
the case m > 2, the claims both on the "symbols" = F„ of the approximating aggregates (jl.4p 
and on the "symbol" of the limiting functional should be more delicate. In the paper |TlJ, the 
uniform (w.r.t. n) analogue of the following "dimensional" condition on the symbol was used: 

(1.10) 3CGM+,7>m-2: fi{{y\\\y - x\\ < r}) < Cr^ , r > 0. 

For the Brownian motion in M™, the following criterium is well known (^, Chapter 8): measure /x 
is a VF-measure iff 

(1.11) sup / w{\\y - x\\)fj.{dy) < +00, 

xGM™ J\\ij-x\\<l 

{max(— Inr, 1), m = 2 
r^-"^, m > 2. 

It is easy to verify that the condition (jl.lOp is sufficient for the condition (jl.lip to hold true. However, 
it is not a necessary one (see Example 15. ip . This, in particular, means that in the main limit theorem 
of [Tl] only the functionals from some proper subclass of the class of VF-functionals (namely, the 
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functionals with their VF-measures satisfying the "dimensional" condition p.lOp ) can be obtained as 
a Umiting ones. 

The main result of the present paper (Theorem 12. ip establishes the weak convergence of the 
functionals (jl.4p for difference approximations Xn of multidimensional uniformly non-degenerate 
diffusions X. In our framework, the class of the difference approximations Xn is wide enough. We 
claim the densities of the transition probabilities for Xn to satisfy a proper version of the local limit 
theorem (property B4), Chapter 4 below). We rely on the results of the papers [13], [T3] while giving 
conditions, sufficient for such a claim. 

The condition, imposed in Theorem 12. II on the "symbols" of the approximating aggregates, is the 
uniform (w.r.t. n) analogue of the condition 

(1.12) lim sup / w{\\y - x\\)n{dy) = 0. 

SIO i^gRm J\\y-x\\<S 

Condition p.lUp is sufficient, but not necessary one for the condition (|1.12p to hold true (see Example 
15. ip . Let us discuss the relation between conditions (jl.l2p and (jl.lip in a more details. In the present 
paper, in order to make exposition more short and transparent, we consider the measures /x with a 
compact supports, only. For such a measure, using the standard estimate for the transition density of 
a non-degenerate diffusion (see [2j, Appendix, §6 and references there), one can check, in a standard 
way, that the condition (jl.lip is also the necessary and sufficient condition for a measure to be a 
VF-measure for uniformly non-degenerate diffusion. The condition (|1.12p is clarified by the following 
statement. 

Proposition 1.1. Let X be a diffusion process, valued in M"^,m > 2, with its coefficients satisfying 
condition Al) of Theorem \2.1i given below. Let fi be a W -measure with a compact support, ip be 
the W -functional, corresponding to this measure and f be its characteristics. Then the following 
statements are equivalent: 

(i) for any t G M"*", the function /*(•) is uniformly continuous on M™; 

(ii) sup^-gigm f\x) ^ 0, 5 i 0; 
(a) ^ satisfies the condition 

We prove Proposition 1 1 . 1 1 in the Chapter 4. One can interpret the statement of this Proposition in 
the following way: any VF-measure, satisfying (|1.12p . correspond to the W^-functional that is regular 
w.r.t. the phase variable. The class of functionals, that can be obtained as the limit ones in the 
context of the main result the present paper (Theorem 12. ip . exactly coincides with the class of the 
functionals, regular w.r.t. the phase variable in a sense given by Proposition 11.11 (see Remark 12. 3p . 

2. The main statement 

We consider the objects, given by (jl.ip - (II. 4p . for m>2. We use notation || • || for the Euclidean 
norm, not indicating explicitly the space this norm is written for. The classes of functions, that 
have k continuous derivatives, and functions, that are continuous and bounded together with their 
k derivatives, are denoted by and C^, correspondingly. The derivative (the gradient) is denoted 
by V. The weak convergence of the (not necessary probability) measures //^ to ji means, by the 
definition, convergence f dfin — > f^m f dfin for every / G C;,(]R'"). 
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Let us formulate the main conditions on the objects involved to (II. 1|) - ()1.4p . 

Al) a G Cfe2^M™,]R"),6 G C^^^jgm^ j^mxm) g^igi- positive constants c, C such that 

c\\ef < {b{x)b*{x)e,e)Mrn < c\\ef, x,ee r"". 

Furthermore, the function satisfies the Holder condition with some positive exponent. 
A2) I.i.d. random vectors {(,k} Eire centered and have the identity for covariance matrix. 
A3) Random vectors {£,k} possess the distribution density p G C^(M™). There exists a function 
: ^ M+ such that 

sup ^(x) < +00, / c^a; < +00 

and 

\\[V'p]{x)\\ <iPix), xeR"", i = 0,...,4. 
A4) Fn{x) > 0,2; G M*",?! > 1 and i sup ^ 0,n ^ 00. 

A5) Measures ^n{dx) = Fn{x)X^{dx) weakly converge to the finite measure /U, that has a compact 
support. 

A6) The following uniform analogue of the condition (11.121) holds true: 

limlimsup sup / w{\\y — x\\)fin{dy) —>■ 0. 

SIO x&M."^ J\\y-x\\<S 

Remark 2.1. If conditions A5) and A6) hold true, then the measure /i satisfies condition (|1.12|) . In 
particular, // is a VF-measure. 

Remark 2.2. The function w{-) is bounded on [5, +00) for any given 6 > 0. Thus, conditions A5) 
and A6) imply that limsup„ ;U„(M™') < +00 and 

limsup sup / w {\\y — x\\)fj,n{dy) < +CO. 

Let us proceed with the formulation of the main statement. Together with the functionals ipn, that 
are piece-wise constant w.r.t. both time variables, we consider the random broken lines, constructed 
from these functionals: 

^n* = + {ns -j + 1)93^'" +{nt-k + l)(p~'", s G 

We interpret the random broken lines ipn as the random elements, taking values in C(T,M^), where 
T = {(s, t)|0 < s < t}. The T4^-functional = ip{X) of the process X we define by the formula (jl.Sp 
(the measure is taken from the condition A5)). 

Theorem 2.1. Let the conditions Al) - A6) hold. Then (X„,'0„(X„)) ^ {X,ip{X)) in a sense of 
weak convergence in C(]R+,M™) x C(T,M+). 

We prove Theorem 12. II in the Chapter 4. 
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Remark 2.3. Let an arbitrary ly-measure fi, satisfying condition ()1.12p . be given. Then one can 
construct a sequence of the functions {-F™} in such a way that conditions A4)-A6) hold true. For 
instance, one can define -F„ by 

Fr,{x)=n - 1 qt{x,y)n{dy)dt, x^W, 

Jo JR"^ 

where qt{x, y) = (27rt)^'2' exp ^— ^^^^ ^ is the transition probabiUty density for the Brownian motion 
W . Then is equal to the value at the point ^ of the characteristics of the Ty-functional of the 
process W , 

and the properties A4)-A6) can be proved analogously to the proof of Proposition 11.11 (see Chapter 
4). 

3. Weak convergence of additive functionals of a sequence of Markov chains 

Our proof of Theorem 12.11 is based on the general theorem on convergence in distribution of a 
sequence of additive functionals of Markov chains, given in the paper [I5j. In this chapter, we give 
a detailed exposition of the objects and auxiliary notions, that are used in this theorem. 

In this chapter we suppose that the processes Xn{-), X{-) are defined and take their values in 
a locally compact metric space (X, p). We say that the process X possesses the Markov property at 
the time moment s G w.r.t. filtration {Qt,t £ if X is adapted with this filtration and for 

every A; G N, ti, . . . , > s there exists a probability kernel {Psti...t^.{x, A), x £ X,A £ B{X'')} such 
that 

(3.1) E[lA{{X{h),...,X{tkmgs] = Pst,...tjX{s),A) a.s., Ael3{X''). 

The measure Psti...tk{x, ■) has a natural interpretation as the conditional finite-dimensional dis- 
tribution of X at the points ti,...,tk under condition {X{s) = x}; below, we use notation 
Pst,...t,{x,-)=P{{X(ti),...,X{tk)) G ■\X{s) = x). 

Everywhere below we claim the process X to possess the Markov property w.r.t. its canonical 
filtration at every s G M"*", and every processes X„ to possess this property (w.r.t. its canonical 
filtrations) at the points of the type ^, i G Z+; this means that every process X„ is, in fact, a Markov 
chain with the time scale, proportional to ^. 

Let the additive functionals ipn be given by the formula (|1.4|) . For the functional <pn, its charac- 
teristics /„, (the analogue of the characteristics of a VF-functional) is defined by the formula 

(3.2) f^'\x)= E[ipf^\Xn)\Xn{s) = x], s = -,ieZ+,t> s,xeX. 

n 

Let us note that the process Xn possesses the Markov property w.r.t. its canonical filtration at the 
points s = ^ , i G Z_|_ and the functional (|1.4|) is the function of the values of X„ at the finite family 
of teh time moments. Therefore the mean value in (j3.2p is well defined as the integral w.r.t. family 
of the conditional finite-dimensional distributions {Pst^...tj,(x, ■),ti, . . . ,tk > s, k £ 'N} of the process 



DIFFERENCE APPROXIMATION FOR LOCAL TIMES OF MULTIDIMENSIONAL DIFFUSIONS 



7 



The fohowing result ([15j. Theorem 1) is an analogue of the well known theorem by E.B.Dynkin, 
that describes convergence of I^-functionals in the terms of their characteristics ([2], Theorem 6.4). 

Theorem 3.1. Let the sequence of the processes Xn be given, providing Markov approximation for 
the homogeneous Markov process X (see Definition \3.1\ below), and let the sequence {(pn = fn{Xn)} 
be defined by Suppose that the following conditions hold true: 

(1) The functions ^Fn{-) are non-negative, bounded on X and uniformly converge to zero: 

6{Fn) *== — supF„(x) — > 0, n ^ oo. 

(2) There exists a function f, that is a characteristics of a certain W -functional (p = 'p>{X) of 
the limiting process X , such that, for every T € M"*", 

sup sup \f^\x) - f*~'{x)\ ^0, n ^ oo. 

s=±,t&{s,T) a:^eX 

(3) The limiting function f is continuous w.r.t. variable x, locally uniformly w.r.t. time variable, 
i.e., for every T G M"^, 

sup|/*(x)-/*(y)| ^0, ||x-y||^0. 
t<T 

Then, for the random broken lines ipn, corresponding to ipn, 

MXn) ^{X) = {ip''\X),{s,t) G T} 

in a sense of weak convergence in C(T,M+). 

//, additionally, Xn =^ X C(M+,X), then {Xn,ipniXn)) =^ {X,ip{X)) n a sense of weak conver- 
gence in C(M+,X) X C(T,M+). 

The notion of Markov approximation, introduced in |l6j , is a key one in Theorem 13.11 

Definition 3.1. The sequence of the processes provides the Markov approximation for the 

Markov process X, if for every 7 > 0, T < +00 there exist a constant K{'y, T) G N and a sequence of 
two-component processes {Yn = {Xn,X'^)}, possibly defined on another probability space, such that 

(i) X„ = X„,X-=X; 

(ii) the processes y„,X„,X" possess the Markov property at the points G N w.r.t. the 
filtration {J^^ = a{Yn{s), s < t)}; 

(iii) limsup_^^P (sup^<^ p (l„ (^) [^)) > 7) < 7- 

The following result, on the one hand, provides an example, that clarifies the given above definition, 
and, on the other hand, gives the opportunity to apply Theorem 13.11 in order to prove the main 
statement of the paper. 

Lemma 3.1. Let a, b be bounded and satisfy Lipschitz condition, m > 1 be arbitrary and i.i.d. random 
vectors {^k} centered, have the identity for covariance matrix and satisfy condition < 
-|-oo for some (5 > 0. 

Then the sequence of the processes Xn, given by provides the Markov approximation 

for the process X, given by the equation (QTTp. 
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At the Example 3 of the paper [15j, the statement of the Lemma was proved with the use of the 
pathwise uniqueness property of the equation (II. ip . The reasonings of such a kind are a quahtative 
ones, and can not provide exphcit estimates for the rate of convergence. Therefore, here we give 
another straightforward proof, that gives possibihty for the further estimates and generahzations. 

Proof. We start from the construction, described in the proof of Theorem 1 [16]. Denote Sn = 
EfeLiCfc! due to CLT, n'^Sn =^ W{1). Condition < +oo ensures that the family {^} is 

uniformly integrable, and therefore the Wasserstein distance between the distributions of the vectors 
N~2Sn' and tends to as — > oo. This means that, for any e > 0, there exist A^ G N and 

random vector (%,Ce) such that 

df Sn, 



/a: 



Ce = W{l). 



Let e > be fixed; we construct the probability space {Q^,J-^,P^) in the following way: 17^ 
(M™)^- X C([0, 1]), = B{n^). Denote the coordinates of a point uj^ £ hy x = ixi, ■ ■ ■ ,Xn,) 
(R™)^=,99 G C([0, 1]). Define the following measures: Q{du,dv) is the joint distribution of (??e,Ce 



Us{dx,u) is the conditional distribution of {.^i, . . . ,^nA under condition {- 

V 

is the conditional distribution of W{-) under condition {1^(1) = v}. We put 



u}, and Vs{d(p, v) 



P\A) 



.)2 



Ue{dx,u)V^{dLp,v) 



Q{du,dv), AeJ'K 



Now we define the probability space {i^,J-,P) as the infinite product of the copies of (17^, JT^, P^). 
For uj = (x^, ip^ ,x^ ,ip^ , . . .) G define the sequence {Ca;(u;), > 1} by the formula 

and the process t G M+} by the formula 

[nt] 



W-{t){co) = -= 



nt]+l 



t - [nt] 



k=l 



n 



t > 0. 



By the construction, the sequence {^fc} has the same distribution with the sequence {^k} and the 
process is a Brownian motion in M™. Now, let us define processes A„,A" by the formulae 
p.2|) . p.3|) and (|l.ip . with {S,k} replaced by {^fc} and W replaced by W""; by the construction, the 
process Yn = (A^) A") satisfies the condition (i) of Definition 13.11 
Also, by the construction, the sets 



■=1 



IN, 



n 



■IN. (i+l)iVei f ' ^ 0, 1, . . . 



are mutually independent. According to ()1.2p . ()1.3p and (II. ip . the value of the process at the given 
time moment i G N, is a functional of H^, . . . , and the values of Y^ A„ or A'^, at any time 
moment t > are a functionals of H^, Hj+i, ... and F„ (^) , A„ (^) or A" (^) , respectively 
Thus, the processes y„,A„ and A" possess the Markov property w.r.t. filtration {jFt^^t G 
generated by y„, at the time moments G N. 
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Let us proceed with the estimation of the distance between X„ and X^. In order to shorten expo- 
sition, we will give the estimate in the partial case m = 1, a = 0; in general case, the argumentation 
is completely analogous, but the calculations take more place. Denote 



1=1 



n 



n 



A," = X^{U) - = / b [X^{s)) dW'is), i e N. 

By the construction, i?[Aj_„|.F(._ J = £'[A"|jr(^_J = 0, therefore, 

(3.3) £;(X„(t,)-^"(t.))' = i^(^n(t^-i)-^"(t^-i))' + i^(Ai,„- A^', ieN. 

We write the decomposition 



A,: 



l)A'' I 



i,n 



E 

1=1 



l-l 



n 



b{Xn{ti-l) 



^{i-l)N,+l 



fn 



The pair of the processes G e N}) has the same distribution with the pair 

G M+},{^fc,/c G N}). Thus, for every A; > 1, the random variable does not depend on 
the values of the process X„ on the interval [0, ^^]. Using this, and taking into account that the 
function h is bounded by a constant B and satisfies the Lipschitz condition with a constant L, we 
get the estimate 



L 



2 



(3.4) ETl^ ^ - E ^ 



1=1 



Xn + 



/ - 1 



n 



Xn (ti-l) 



2 ^2N^ 



< — yB^ < ^ 



1=1 



2n^ 



Let us also write the decomposition 

A^ = b (X" (t,_i)) • [W^^\U) - W^"(t._i)" 
where the second summand (we do not write it explicitly) can be estimated analogously to ()3.4p : 



(3.5) 

At last, we write the decomposition 



r2R2 Ar2 



^ii-l)N^+l 



1=1 ^ 



blX^iU-i)) ■ #"(t,) - H^"(ii-i) 



b(Xn{U-l)) -b(x^{ti^i)) ■ #"(ti) - t^"(t*-i) 



+ \l^-H^n{u^i))-ri„ 
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where = Y^i^=i i{i-i)N,+i - [U - k^i] ^ W'"'iti) - . By the construction, for every 

i > 1, the random variables W"'{ti) — l^"(tj„i) and do not depend on the values of the process 
Yn on the interval [0, 77^]- Moreover, £'[T^]^ < e. Therefore, the following estimates hold true 



(3.6) E blXn - 6 X" W^iU) - 



2 r 



< 



n 



E{Xn{ti-l)-X^{ti.i)f 



(3.7) 



E 



b X„. it 



-L ,f- 



Now, using the decomposition p.3p . the Cauchy inequality (X]fe<4 ^fc)^ — ^ X]a:<4 ^fc ^^"^ estimates 
(I33])-(132I), we obtain 



(3.8) E{Xn{k)-X''{U)f < + 
Iterating (13. Sp . we obtain 

(3.9) ^(X„(t,)-^"(i*))' < 



EiXniU^l) - X'^(t._l))2 + 



+ 



n 



AL'^B'^N, 



n 



n 



Let ti < T, that means that i < jp^. Then the sum in the right hand side of (13. 9p contains at most 



summands, and every summand is not greater than exp 

AL'^B'^N.T 



estimate 



. This provides the 



(3.10) 



E{Xr,{U)-X^{U)Y < 



+AB'^Te 



n 



The sequence {(X„(tj) — X"(tj)), i G N}, by the construction, is a martingale, thus, using the maximal 
martingale inequality ([17], Ch. VII, §3), we obtain the estimate 

^AL'^B^N^T 



(3.11) 



P 



max(X„(t,)-^"(t*))' >7' 

i:t,<T 



<7 



-2 



+ AB^Ts 



n 



e^^'^, 7>0. 



Now we can complete the proof of the Lemma. For a given 7, T choose e > in such a way that 

16i?2^.e^^'^•e<7^ 

and proceed with the construction, described above, with this £. We have already seen that, under 
this construction, conditions (i) and (ii) hold true with K{"f^T) = Nf,. The estimate (|3.1ip provides 
that the condition (iii) holds true with the same K{^,T). The lemma is proved. 

Remark 3.1. Denote by K(7,T) the minimum of the set of such numbers K gN, that here exists a 
process Yn, satisfying conditions (i) - (iii) of Definition 13.11 with K{'y,T) = K. In the paper [16], it 
is shown (the part II of Theorem 1) that, in the basic case a = 0,b = /Mm, 

supK(7, T) = +00 

7>0 

as soon as the distribution of ^1 differs from the normal one. One can say that, while the accuracy 
of the approximation becomes better (the accuracy is described by the parameter 7), the Markov 
properties of the two-component process necessarily become worse (these properties are described by 
K(7,T)). 
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4. The proofs of Theorem 12.11 and Proposition 11.11 

We reduce the proof of Theorem 12.11 to the verification of the conditions of Theorem 13.11 The 
sequence Xn provides the Markov approximation for the process X due to Lemma |3. 11 Condition 1 
of Theorem 13. II holds true due to condition A4). Let us check that the conditions 2 and 3 hold true. 
The characteristics /* of the functional c/^*'* = {X{r)) dr has the form 

f\x) = [ [ Prix, y)fi{dy) dr, x G M"^, t > 0, 
Jo JR™- 

where {pt{x,y),t > 0, x,y £ M™} is the transition probability density for the process X. Existence 
of such a density under condition Al) is a standard result of the theory of parabolic equations. 
Moreover, this density possesses the following properties (see, for instance. Appendix to [2], §6, and 
references there). 

Bl) The function p : (t,x,y) i— > pt{x,y) is uniformly continuous on [S,T] x x M"* for every 
0<6 <T. 

B2) There exist a constants M, a > such that 

, , / ally — a^P' 
pt{x, y) <Mt 2 exp ( ^ 

B3) There exist a constants Mi, M2, a, /?, A > such that 

,^ ™ / ally — a;|P\ ,^ m,), / /5||w — x|P 
pt{x, y) > Mit-^ exp f ^ " j - exp f- ^"^^ " 

Under conditions A2),A3), the processes X^ possess the transition probability densities p^{x,y) 
at the time moments t G ^N, that means that 

P{Xn{t) (^T\Xn{s) =x) = I p1_Ax,y)dy, s,te-Z,s<t, x G M", F G «B(M™), 

Jr n 

and, moreover, for every t G ^N, the function is a continuous one. The characteristics /„ of the 
functionals ipn can be expressed through these densities by the formula 

(4.1) f/[x) = -Fn{x) + - I F^{y)pl{x,y)dy = f/-'{x), s<t, 

Below, we denote /n'* = Theorem 2.1 [13] and Theorem 1 [TJj imply that, under conditions Al) 
- A3) , the following estimate for the deviation of the densities p^ from the limiting density pt holds 
true. 

B4) For any T G M+, 

sup sup ^/n-t~^ ■ (1 + ( ^^^ ] ] ■\pt{x,y) - p^{x,y)\ < +00. 
Denote, for (5 > 0, 

fl{x) = fix) - f^'ix) = f I Prix, y)piidy) dr, x^W^,t>Q, 



X G R'". 
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f^,{x) = fi{x)-C\x) = - I Fn{y)pl{x,y)dy, s<t,xGR^. 

Conditions A4), A5) and properties B1),B2),B4) imply the following statement. 

Proposition 4.1. For any 6 > 0, 

(i) the function fs : {t,x) i-^ fsi^) uniformly continuous on [0,T] x M™; 

(ii) the functions fs^n ■ {t,x) i-^ fns(^) converge, as n ^ oo, to the function fs uniformly on 
[0,r]xM"^. 

The proof of Proposition 14. II is quite standard and we omit it here (see, for instance, the arguments 
given in the proof of Theorem 3 [E]). It follows from Proposition 14. 1 1 that . in order to prove that the 
conditions 2 and 3 of Theorem 13.11 hold true (and therefore, to prove the Theorem 12. 1|) . it is sufficient 
to prove the following two relations: 



(4.2) 



lim sup / (x) = 0, 



(4.3) 



lim lim sup sup /„(x) = 0. 



Let us prove ()4.3|) . In the exposition below, we suppose that 6 < 1. It follows from B2),B4) that, for 
any T G there exists a constant Ct £ such that 



pnx,y)<CT-t-^ 



exp 



a\\y — x\\ 



+ 1 + 



\y - x\ 



m\ — 1 



Then the formula (j4.1|) implies the estimate 

(4.4) sup fi{x) < - sup Fn{x)+ sup / Kn^siWy - x\\)^JLn{dy) 



i e -N,t < r,n G N. 
n 



where 



k=l 



n 



exp -az •t +1+ -^'a/t 



m\ —V 



z G 



For any A; G N and any i G |^^, , the following inequalities hold true: 



2 I k + l\~ 
< I — — ) -t-- < 2- -t--, 



m m m 



n 



exp ( -az ■ -j^] + \'^ + \ z ■ \ j — 



m\ —1 



< 



exp 



az^\ / f z 



m\ —1 



and therefore, for any k gN, 
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The estimate (I4.5p implies the following estimate for the kernel Kn^s'- 



71 

< Ks{z) = 2fCi ■ / t- 



exp 



az 



+ 1 + 



z 

7t 



m\ — 1 



dt < 



exp 



az 



+ 1 + 



z 

Vt 



m\ — 1 



dt 



We estimate the kernel using the calculations, analogous to those made in [2j, §6.1. Making the 

2 

change of variables u = 4-, we get the formula 



(4.6) Ksiz) = 2^C\-z^-"'-G(^ 



G{y) 



U 2 



ex.p{—au) + I 1 + u 2 



du. 



The following properties of the function G can be verified straightforwardly: 

G(-) is bounded, m > 2 



(4.7) 



G{y) ^0, +0O, 



G(y) ~ -21ny, y ^0+, m = 2 



Let us give two estimates. First, 
(4.8) 



Ksiz) < D ■ w{z) ■ G{5- 



z > 6* 



(here and below we denote by D a positive constant such that its explicit value in not needed for us; 
the concrete values of D may vary from line to line). For m > 2, this estimate follows immediately 
form the representation (j4.6p : for m = 2, one should recall additionally that (for m = 2) w{z) > 1. 
Second, 



(4.9) 



Ks{z) < Dw{z) 



z G 



For m = 3 this estimate holds true since the function G is bounded, for m = 2 it holds true due to 
the estimates 



Giy)<D-w{y), G 



z^\ „ /z^ 



< D-w 



< D ■ w(z^) < 2D ■ w(z) 



5 J - \6 

(we used here that 5 < 1). The formula (14. 4p . condition A4) and estimates (I4.8p . ()4.9p imply 
limsup sup fn{x) < D • G(5~2) . limsup sup / w{\\y — x\\)fin{dy)+ 

n—>+oo leK™ n^+oo xSK™ Jr™ 



(4.10) 



+ D • limsup sup / ^ w{\\y - x\\)iin{dy). 

n^+oo zeR™ J\\y-x\\<S^ 



Both summands in the right hand side of the inequality (I4.10p tend to zero as (5 | 0. This follows from 
(14. 7|) . Remark 12.21 and condition A6). Thus, the relation (14. 3p is proved. The proof of the relation 
(|4.2p is completely analogous and we omit it here. 

We have verified that the processes Xn,X and the functionals defined by (jl.ip - ()1.4p . under 
conditions Al) - A6), satisfy all the conditions of Theorem 13.11 Using this theorem, we obtain the 
statement of Theorem 12.11 
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The proof of Proposition f777l The family {Z*^} is monotonous w.r.t. 6. Therefore, the imphcation 
(i)=^(ii) fohows form the Dini's Theorem and the fact that, due to compactness of the support of ^ 
and the property B2), for any given t G M+, 

f^{x) 0, ||x|| +00. 

The uniform hmit of a uniformly continuous functions is a uniformly continuous function; this, to- 
gether with Proposition 14. 11 provides the implication (ii)=^>(i). The implication (iii)=^>(ii) is contained 
in the relation (j4.2p . In order to prove the inverse implication, we use the property B3) and calcula- 
tions, analogous to those made above. We write 



f'{x)> I RsiWy - x\\)i,{dy), Rs{z) ^ 



Mir~ exp M2t"~+'^exp - — 







t \ t 



dt. 



2 



Making the change of the variables w = 4-, we get 



'2 



6 



poo poo 

Gi{y) = Ml u^~'^exp[-au]du, G2(y) = M2 / u^'"^'^ exp[-Pu]du. 

Without losing generality, we can suppose that A < 1. Let us consider two cases. If m > 3, then the 
functions Gi,G2 are bounded, and, for z <5, 

Rsiz) > z^-™ • [Gi{5) - 6'^^G2{S)] > w{z) ■ D 

with 6 small enough. Thus, for small enough 6, 

(4.11) f\x)>Df w{\\y-x\\)ii{dy). 

J \\y—x\\<5 

If m = 2, then Gi{y) ~ — Iny, y J, 0, and therefore, for small enough 5 and z < S, 

Gi(^^ > Dlnir > Dln- = Dwiz). 
\ J z 

Next, G2{y) ^ \y^^^ y i 0, and thus, for small enough 5 and z < 6, 



5 ) - \z^ 

The two latter estimates and the fact that w{-) > 1 for m = 2 provide that the estimate (14. lip holds 
true in the case m = 2, also. This estimate proves the implication (ii)=^(iii). The proposition is 
proved. 

5. Examples 

In order to illustrate the domain of application of Theorem 12. H we give two examples. In our 
first example, we construct the measure fi, satisfying condition (|1.12p . but not satisfying (|l.lUp . The 
ly-functional, that corresponds to this measure, can not occur as a limiting one in the framework 
of the paper [11]; on the other hand, this functional belongs to the class of the limiting functionals, 
obtained via Theorem 12.11 (see Remark 12. Sp . 
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Example 5.1. Let m = 2, {r^, A; > 1} be a sequence of positive numbers (it will be defined precisely 
later on), and the measure fi to have the form fi = ^^yiQk ■ fJ-k, where {Qk,k > 1} is a certain 
weight sequence and fi^ = ^Sk is the surface measure on the circle Sk = {?/|||y|| = ri^},k > 1 (the 
measures are normalized in such a way that //fe(S'fc) = l,k > 1). 
For any x G M^, 5 € (0, e~^), denote B{x, 6) = {y| \\y — x\\ < 6}, then 

/ wi\\y-x\\Mdy) = J2Qk-Viix), F/(-) = /" wi\\y - ■]]) fik{dy). 

A\y-^\\<^ k JSknB{x,5) 

Note that, up to a term is the simple layer potential, generated by the measure lB{x,S)dfJ'k: 

concentrated on Sk (see, for instance, fl8], §21). This potential is dominated by the potential Vk 
generated by the measure ^k, and both these potentials are a continuous functions, harmonic in both 
{y|||y|| < fk} and {y|||y|| > r^}. The maximum principle provides that 
(5.1) 

Vk{z) = Vk{0) = In (J-^ , ||z|| < rfc, V^ix) < Vk{x) < In (J-^ , x G M^ 5 G (0, e~'), k > 1. 

Due to the maximum principle, takes its maximum value on the circle Sj.. If 6 < r^, then 
Sk n B{x, 5) is the arch on the circle Sk, and it is easy to verify that the corresponding maximal value 
is taken at the middle of this arch. This reasoning and the straightforward calculation, that is easy 
and omitted, give the estimate 

(5.2) v,\x)<D-(^)-ln^^ 



Now, we put rfc = 2^'^^, = k~'^, k > 1. Let 5 G (0, e^^) be fixed and x G M™" be arbitrary, let us 
estimate J\^y_x\\<s^i\\y ~ ^IDf^idy)- For < 25, we have B{x,6) C B{0,36), an therefore, due to 

(EH), 

(5.3) I my - .iDM^i.) < E (;^) = E ^ s " ■ [i" (s)] " ■ 

If ||x|| > 2(5, then there exists at most one value k = k^, such that B{x, 5) n Sk^ 7^ 0, it being known 
that rfc. > 5. Then, due to (pT]) . (j5:2]) . 

(5.4) / w{\\y-x\\)^i{dy) < D- sup [\nrk\~^-(-\ -In (-\ = D- sup L (-\\ , 

J\\y-x\\<S k:rk>S K'TkJ \0 J k:rk>S ^{Tk) V \0 J _ 

where the function <I>(t) = t[lnt]^,t > is non-decreasing in some neighborhood of zero. From 
(1E3D,(153D, the estimate 



sup / w{\\y - x\\)n{dy) < D 

xeK™ J\\y-x\\<5 

follows, that provides (I1.12p . On the other hand, 

fi(B(0,ri))>(3t = 



mil 



In Tfc 
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and thus, for every 7 > 0, 



limsupr n{B{0,r)) = +00, 

rlO 



that means that condition (jl.lOp does not hold true. 

In the second example, we construct a VF-measure ^, that does not satisfy condition (11.120 : the 
M^-functional, corresponding to this measure, can not be obtained in Theorem 12.11 as a limiting one. 

Example 5.2. Let m = 2, we put = 2-^\ak = (^,0) E M™, Sk = {y\\\y - ak\\ = rk},Qk = 
k^'^, k > 1, flk = surface measure on the circle Sk, normalized in such a way that fj.k{Sk) = 

1,A: > 1. We put 

k 

and show that /i is a VF-measure that does not satisfy ()1.12p . By the construction, there exists iV G N 
such that, for any x G M™", the relation 



(5.5) ||x — ak\\ < 2 



-Vk 



holds true for at most N values G N. Furthermore, if A; > 4 and (j5.5|) holds true, then ||?/ — x|| < 
e~^,y G Sk- Therefore, for a given x G M™, for such k > A that (j5.5p holds true, using the maximum 
principle we obtain the estimate 

Qk ■ I w{\\y - x\\)flk{dy) < Qk ■ ( In {- — ^- — -A fikidy) = Qk ■ hi ( — ) =Qk-k'^ = '^- 



\\y-x\\<l JSk \\\y~^k\\/ \rk 

If k > 1 and ()5.5p fails, then \\y — x\\ > 2^^^^,y G Sk- Thus, for a given x G M"*, for k > 2 such 
that (15. 5p fails, we have the estimate 

Qk ■ I w{\\y - x\\)fLk{dy) < Qk ■ {Vk + 1) = /c"! + k-\ 

Furthermore, every measure flk is a M^-measure, and thus 

(Tfe = sup / w{\\y - x\\)ilk{dy) < +00, k>l. 

xeM™ J\\y~x\\<i 

The three latter estimates imply that 

r ^ 3 

sup / w{\\y -x\\)li{dy)<y2k-^ak + N + y2{k- 2 +k-^)< +00, 

xm^J\\y-x\\<l fc^^ ^,>5 

that means that /i is a T4^-measure. On the other hand, for 5k = rk,Xk = a-k we have that 

/ w{\\y-Xk\\)il{dy)>Qk- j w{\\y - Xk\\)ilk{dy) = I 

and 5k — > 0, A; ^ 00. This means that, for /i, condition p.l2p fails. 
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